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1 Introduction 

It is assumed in the celebrated Black-Scholes model that the volatility of a 
stock is constant. However, empirical studies do not support this assumption. 
In more recent models, the volatility of a stock is represented by a stochastic 
process. Well-known examples of stochastic volatility models are the Hull- 
White, the Stein-Stein, and the Heston model. The volatility processes in 
these models are a geometric Brownian motion, the absolute value of an 
Ornstein-Uhlenbeck process, and a Cox-Ingersoll-Ross process, respectively. 
For more information on stochastic volatility models, see ^ and [B]. 

A stock price model with stochastic volatility is called uncorrelated if 
standard Brownian motions driving the stock price equation and the volatility 
equation are independent. In [7], [9], and [10], sharp asymptotic formulas 
were found for the distribution density of the stock price in uncorrelated 
Hull- White, Stein-Stein, and Heston models. Various applications of these 
formulas were given in [8] and [TI]. The results obtained in [9] and [10] will 
be used in the present paper. 

It is known that the stock price distribution density in an uncorrelated 
stochastic volatility model possesses a certain structural symmetry (see for- 
mula flT^ below). This implies a similar symmetry in the Black-Scholes 
implied volatility, which does not explain the volatility skew observed in 
practice. To improve the performance of an uncorrelated model, one can 
either assume that the stock price process and the volatility process are cor- 
related, or add a jump component to the stock price equation or to the 
volatility equation. The stock price distribution in the resulting model fits 
the empirical stock price distribution better than in the uncorrelated case. 
However, passing to a correlated model or adding a jump component may 
sometimes lead to similar effects or may have different consequences (see e.g. 
m and [2j). Examples of stock price models with jumps can be found in [3], 
|15j . and [TB]. We refer the reader to [4] for more information about stock 
price models with jumps. An interesting discussion of the effect of adding 
jumps to the Heston model in contained in [Ti] . 

An important jump-diffusion model was introduced and studied by Kou 
(see [15j and [I6]). This model can be described as a perturbation of the 
Black-Scholes model by a compound Poisson process with double-exponential 



law for the jump amplitudes. In the present paper, we consider similar 
perturbations of stochastic volatility models. Our main goal is to determine 
whether significant changes may occur in the tail behavior of the stock price 
distribution after such a perturbation. We show that the answer depends 
on the relations between the parameters defining the original model and the 
characteristics of the jump process. For instance, no significant changes occur 
in the behavior of the distribution density of the stock price in a perturbed 
Heston or Stein-Stein model if the value of the parameter characterizing the 
right tail of the double exponential law is large. On the other hand, if this 
value is small, then the distribution density of the stock price in the perturbed 
model decreases slower than in the original model. For the Hull- White model, 
there are no significant changes in the tail behavior of the stock price density, 
since this density decays extremely slowly. 

We will next briefly overview the structure of the present paper. In Sec- 
tion |2l we describe classical stochastic volatility models and their perturba- 
tions by a compound Poisson process. In Section [3] we formulate the main 
results of the paper and discuss what follows from them. Finally, in Section 
m we prove the theorems formulated in Section [31 

2 Preliminaries 

In the present paper, we consider perturbations of uncorrelated Stein-Stein, 
Heston, and Hull- White models by compound Poisson processes. Our goal 
is to determine whether the behavior of the stock price distribution density 
in the original models changes after such a perturbation. 

The stock price process X and the volatility process Y in the Stein-Stein 
model satisfy the following system of stochastic differential equations: 

dXt = ^jiXtdt+\Y,\X,dWt , 

dYt = q{m-Yt)dt + adZt. ^^ 

This model was introduced and studied in [18]. The process Y, solving the 
second equation in ([1]), is called an Ornstein-Uhlenbeck process. We assume 
that // G M, g > 0, m > 0, and a > 0. 

The Heston model was developed in [12]. In this model, the processes X 
and Y satisfy 

' dXt = fiXtdt + ^/YtXtdWt 
dYt = q{m-Yt)dt + c^tdZt , ^ ' 



where /z G M, g > 0, m > 0, and c > 0. The volatihty equation in ([2]) is 
uniquely solvable in the strong sense, and the solution y is a non-negative 
stochastic process. This process is called a Cox-IngersoU-Ross process. 

The stock price process X and the volatility process Y in the Hull- White 
model are determined from the following system of stochastic differential 
equations: 

' dXt = fiXtdt + YtXtdWt . s 

dVt = uYtdt + CYtdZt. ^^ 

In ([2]), /i € M, z^ G M, and ^ > 0. The Hull- White model was introduced in 
[13j . The volatility process in this model is a geometric Brownian motion. 

It will be assumed throughout the paper that standard Brownian motions 
W and Z in ([1]), (|2]), and (|3]) are independent. The initial conditions for the 
processes X and Y will be denoted by xq and yo, respectively. 

We will next discuss perturbations of the models defined above by a 
compound Poisson process with jump amplitudes distributed according to a 
double exponential law. Perturbations of the Black-Scholes model by such 
jump processes were studied by Kou in [15] and by Kou and Wang in [T6] . 
Some of the methods developed in ^5\ will be used in the present paper. 

Let A^ be a standard Poisson process with intensity A > 0, and consider 
a compound Poisson process defined by 

Jt = J2{V,-l), t>0, (4) 

where Vi are positive independent identically distributed random variables 
that are independent of {A^t}i>o- ^^ ^^ ^^^'^ assumed that the distribution 
density f oi Ui = log Vi is double exponential, that is, 

f{u) = pr/ie-''^"l{„>o} + gr]2e''^"l{„<o}. (5) 

where 771 > 1, 772 > 0, and p and q are positive numbers such that p + q = 1. 
Consider the following jump-diffusion stochastic volatility models: 

dXt = fiXt- dt+\Yt\ Xt^ dWt + Xt- dJt ,^. 

dYt = q{m- Yt) dt + adZt ^ ^ 

(the perturbed Stein-Stein model), 

dXt = fiXt- dt + VYtXt- dWt + Xt-dJt ,„. 

dYt = q{m- Yt) dt + Cy/YtdZt, ^ ' 



(the perturbed Heston model), and 



dXt = fiXt_dt + YtXt^dWt + Xt^dJt 
dYt = uYtdt + CYtdZt, 



(8) 



(the perturbed Hull- White model). It is assumed in ([6]), ([7]), and ([8]) that the 
compound Poisson process J is independent of standard Brownian motions 
W and Z. 

We will next formulate several results of Gulisashvili and Stein. For the 
uncorrelated Heston model, there exist constants Ai > 0, A2 > 0, and A3 > 2 
such that 



Dt{x) = Ai(logx) 4 



3 1 qm 



c^ e 



'^a;""^' 



8 (l + o({\ogx)-^ 



(9) 



as a: — )■ cxD (see [ID]). For the uncorrelated Stein-Stein model, there exist 
constants Bi > 0, B2 > 0, and B^ > 2 such that 



Dt{x) = Bi(logx)-h^^^^^x-^^ (1 + (aogx)--A] 



(10) 



as X — )■ 00 (see 53). Finally, in the case of the uncorrelated Hull- White 
model, there exist constants 61 > 0, 62 and 63 such that following formula 
holds (see [9] and also Theorem 4.1 in [lOj): 



Dt{x) = bix~ (log x) 2 (log log X 
1 



exp 



log 



2t^2 

1 + ( (log log x)" 2 



1 / 2 log X 

yo 



t 



,63 



+ - log log 



1 / 2 log a; 

1/0 



t 



(11) 



as a; — )■ 00. The constants in formulas (|9]), (ITO|) . and (TTT]) depend on the 
model parameters. Explicit expressions for these constants can be found in 
[9] and [10]. The constants A3 and S3, appearing in ([9]) and (fTOj) . describe 
the rate of the power-type decay of the stock price distribution density in 
the Heston and the Stein-Stein model, respectively. The explicit formulas for 
these constants are as follows: 



A. 



3 V8C + 
2^ 2Vt 



t 



with C 



2^G' + ^4" 



;i2) 



and 



In f lT2|) and flT3|) . r^ denotes the smallest positive root of the entire function 



z h^ z cosz + s smz. 



Formulas ( 1121) and ( 1131) can be found in [TO] . 

The distribution density density Df in uncorrected stochastic volatility 
models satisfies the following symmetry condition: 

£2£!:)'z,,((-^)^AW, x->o. (14) 

(see Section 2 in [ID]). This condition shows that the asymptotic behavior 
of the stock price distribution density near zero is completely determined by 
its behavior near infinity. 

3 Main results 

The following theorems concern the tail behavior of the stock price distribu- 
tion density in perturbed Stein-Stein, Heston, and Hull- White models: 

Theorem 1 Let e > 0. Then there exist ci > 0, C2 > 0, and Xi > such 
that the following estimates hold for the distribution density Df of the stock 
price Xt in the perturbed Heston model: 

ci 1^ + ^^ I < A(a;) <C2 (^ + ^- I (15) 



for all X > Xi- In US]) , the constant A^ is given by u^) and the constants C2 
and Xi depend on e. 

Theorem 2 Let e > 0. Then there exist cs > 0, C4 > 0, and X2 > such 
that the following estimates hold for the distribution density Dt of the stock 
price Xi in the perturbed Heston model: 

03 {x^'-^ + x"^-^) < Dt{x) < C4 {x^'-"^-' + x"^"^-") (16) 

for all < x < X2. Here the constant A^ is the same as in TheoremU\ and 
the constants C4 and X2 depend on e. 



Theorem 3 Let e > 0. Then there exist C5 > 0, Cg > 0, and X3 > such 
that the following estimates hold for the distribution density Dt of the stock 
price Xt in the perturbed Stein-Stein model: 

for all X > x^. In (11), the constant B3 is given by l[T^) and the constants cq 
and X3 depend on e. 

Theorem 4 Let e > 0. Then there exist c-j > 0, Cg > 0, and X4 > such 
that the following estimates hold for the distribution density Dt of the stock 
price Xt in the perturbed Stein-Stein model: 

c-j {x^'-^ + x''^-^) < Dt[x) < cg (x^^-3-" + x''^-^-') (18) 

for all < X < X4. Here the constant B^ is the same as in Theorem\^ and 
the constants Cg and x^ depend on e. 

We will prove Theorems [1]|1] in Section HI In the remaining part of the 
present section, we compare the tail behavior of the stock price distribution 
density before and after perturbation by a compound Poisson process. 

Let us begin with the Heston model. It follows from Theorem [1] that if 
1 + 771 < As, then 

< Dtix) < — r— -— , X > Xi. 



Therefore, formula (jH]) shows that that if the condition 1 + r/i < A3 holds, 
then the tail of the distribution of the stock price in the perturbed Heston 
model is heavier than in the original model. 

On the other hand, if 1 + r^i > ^3, then Theorem [1] implies the following 
estimate: 

< Dtix) < ^. — 7, x>xi. 



Now formula iQ shows that if 1 + ?7i > ^.3, then there are no significant 
changes in the tail behavior of the distribution density of the stock price 
after perturbation. Similar assertions hold for the Stein-Stein model. This 
can be established using Theorem |3] and formula (ITOj) . 



Next, suppose x -^ 0. Then we can compare the behavior of the dis- 
tribution density of the stock price in unperturbed and perturbed models, 
taking into account Theorem [2], Theorem |U formula ([9]), formula fITOl) . and 
the symmetry condition flT^ . For instance, if 772 < As — 2 in the perturbed 
Heston model, then 

for all X < X2. On the other hand ii r]2 > A3 — 2, then 

d3X^'-^ < Dt{x) < diX^'-'^-' 

for all X < X2. Similar results hold for the Stein-Stein model. 

For the Hull- White model, there are no significant changes in the tail 
behavior of the stock price distribution after perturbation. This statement 
follows from the assumption rji > 1 and from the fact that the stock price 
density in the unperturbed Hull- White model decays like x'"^ (see formula 

m)- 

4 Proofs of the main results 

The proofs of Theorems [Hll] are based on an explicit formula for the distri- 
bution density Df of the stock price Xt in perturbed Heston, Stein-Stein, 
and Hull- White models (see formula (1221) below). Note that the stock price 
process X in the perturbed Stein-Stein and Hull- White models is given by 

Xt = xoexp J /zt - i /" Y^^ds + f YsdWs + Y.uA , (19) 



while for the perturbed for Heston model we have 

Xt = Xoexplfit-- Ysds + / V^dlV, + J2^^ 



(20) 

Formulas flT9|) and (|20|) can be established using the Doleans-Dade formula 
(see, for example, [H]). We will denote by /it the distribution of the random 
variable Jt defined in (j3]). It is not hard to see that the following formula 
holds: 



oo „ 

^t{A) = 7ro(5o(A) + Vtt,, / f<^\u)du (21) 

n=l -^A 

where ttq = e~'^*, 7r„ = e~'*'*(n!)~^(At)" for n > 1, A is a Borel subset of R, 
and / is given by ([5]). The star in fl2T]) denotes the convolution. 

The distribution density D^ of the stock price Xt in uncorrelated models 
of our interest is related to the law of the following random variable: 

for the Stein-Stein and the Hull- White model, and 



^i 



for the Heston model (see [9] and [10]). The distribution density of the 
random variable at is called the mixing distribution density and is denoted 
by rrit. We refer the reader to [H], [ID], and [IB] for more information on the 
mixing distribution density. 

The next lemma establishes a relation between the mixing distribution 
density rrit in the uncorrelated model and the distribution density Dt of the 
stock price Xt in the corresponding perturbed model. 

Lemma 5 The density Dt in perturbed Stein-Stein, Heston and Hull-White 
models is given by the following formula: 



where rrit is the mixing distribution density and fit is defined by l[21\) . 

Proof: We will prove Lemma \5\ for the Heston model. The proof for the 
Stein-Stein and the Hull- White model is similar. For the latter models, we 
use formula (fT9|) instead of formula (120|) . 

Put Tt = X]i=i ^i- Then for any r/ > 0, formula (120|) gives 




F{Xt<rj 



P 



E 



V , ta^ 



y/YJW, + Z < log ^- + 

2, POO -| 



oo ^ — oo 



2iTtai 



exp 



2ta2 



fitidu)dz, 



where z^, = log H -. Making the substitution z = log H -, 



we obtain 



XqC' 



fit 



xoe'^* 



P (Xt < v) 



E 



rj foo 



1 



taf 



oo J -oo V 27rtoj 
exp 



exp 



2ta2 



ixt{du) — 

X 



^ r^ roo I {\og-^ + iM^ - uf [ ^^^^^^^mt{y) ^Jx 



fitidu) 



-dy- 



-ooJO J-oo I 2t|/2 I ■ ^^ ' ^/2'Kty " X 

It is clear that the previous equality implies Lemma [H 
Remark 6 It follows from Lemma [5] that 



Dt{x) 



J 3~ 

V27rtx2 



e^lit[du) exp <^ 'j—- ^(^V- (22) 



y 



2ty^ 



This representation will be used below to obtain two-sided estimates for 
the distribution density of the stock price in perturbed stochastic volatility 
models. 

Proof of Theorem [T], The next lemma will be needed in the proof of 
Theorem [H 
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Lemma 7 Let f be the density of the double exponential law (see formula 
^). Then for every n > 1, the following formula holds: 



/*(")(t,) = e-''^«^P„ 



kVi 



k=l 

n 



(k-l) 



+ e"^" Y. ^". 



kV2 



fc=l 



{k-l)\ 



k—1 n 

■u 1{m>o} 



-u) ]1{«<0}, 



(23) 



where 



n-l 



Pn,k - 2^ 



i=k 



n — k — 1\ (n 
i — k 



i—k / \ n—i 

Vl \ I V2 \ i n-i 



^1+^2 J V??l+??2 



for all 1 < k < n — 1, and 



n-l 



^n,k / ^ 



i=k 



n — k — l\ fn\ f rji 



i-k J \ij \r]i +r]2 



n—i / \ i—k 

V2 \ ^n-i^i 



Vl +V2 



for all 1 < k < n — 1. In addition, Pn^n = P" ond Qn,n = (p ■ 

Lemnia [7] can be established using Proposition B.l in [?] and taking into 
account simple properties of the exponential distribution. 

The next statement follows from Lemma [7] and formula (I^T]) : 

Lemma 8 For every Borel set A cM., 



IJt{A) = TToSoiA) + Gi{u)e~'"''du+ GsHe^^^dw, (24) 

'An[o,oo) J Ar\{-oofi) 



where 



Gi{u) = J2 



fc=0 



k+l oo 



k\ 



S" n P 



k+l 



n=k+l 



U 



(25) 



and 



G2{u) = J2 



fc=0 



, ■ / ^ '^nQn,k+l 
n=k+l 



{-U)\ 



(26) 
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Our next goal is to estimate the rate of growth of the functions Gi and 
G2 defined by (I2S]) and ([2nD. 

Lemma 9 For every e > the function Gi grows slower than the function 
M I— )■ e^" as zx —7- 00. Similarly, the function G2 grows slower than the function 



Proof: We will prove the lemma by comparing the Taylor coefficients 



"''' ^ IT^l^^ X^ TrnPn,k+l, k>0, 



k\ 

n=k+l 

of the function Gi and the Taylor coefficients bk = ^£^, A; > 0, of the function 
e^"". We have ak < hk for k > k^. The previous inequality can be established 
using the estimate 

00 00 

r/l^^ Y^ T^nPn,k+l < Vl^' Yl ''"' 
n=k+l n=k+l 

and taking into account the fast decay of the complementary distribution 
function of the Poisson distribution. 

This completes the proof of Lemma [H] for the function Gi. The proof for 
the function G2 is similar. 
■ 

The following lemma was obtained in [10] (formula (54)): 

Lemma 10 Let rut he the mixing distribution density in the Heston model. 
Then there exist constants Hi > and H2 > 0, depending on the model 
parameters, such that 

°° My) „,.„ f f ^' , ty' 




, ^"Pi"V^ + ~'^^^ 



= i/ic-t+^e^^v^exp {-^^^^] (1 + (o;-^ 

as w —7- 00. The constant C in the previous formula is given by [T^j . 

Proof of the estimate from below in Theorem [Tl We will use 



formula (1221) in the proof. Put z = log ^^t ■ Then we have 

12 



3 / e^ntidu) / exp<^ -— 

^27rtx2 Jr Jo y I 2ty- 



A(2;) = ,-- 3 / e^ntidu) I ^^exp { ^^-^ ^}dy- (27) 



Note that for the uncorrelated Heston model the following formula holds: 

V^^ rm,{y) r z' ty'\ 

Dt{x) = 3 / exp <^ -— ^ - ^ r ^2/ (28) 

v/27ftx2 Jo y I 2ty^ 8 J 

(see [10]). 

Let p be any increasing function of z such that p{z) < z and z — p{z) — )■ oo 
as 2; — 7- oo. Then fl27|) gives 



A(x)>/i + /2, (29) 

where 

h = , - / e2Ht{du) exp<^ — tt C dy (30) 



^27rtx2 Ji Jo y [ '^ty 

and 

e2/ii(du) / exp<^ — ^}dy. (31) 



/27rta;t A ' ^"'^^"""' J^^ y " "^^ (^ 2ty'^ 



Throughout the remaining part of the section, we will denote by a a 
positive constant which may differ from line to line. Since the function Gi is 
increasing on (0, 00) and (12^ and (!25l) hold, we have 

z+l POO ,_, f^\ ( 1 j„,2 



mt{y) ( 1 t?/2 



I2 > ax 2 62 e ''^"d'u / exp < > dy, x > Xq. 

Jz Jo y [ ^W 8 J 

It is known that / y^^mt{y)dy < 00 (see (ID])- Therefore, the second inte- 

Jo 
gral in the previous estimate converges. It follows that 

"Z + l 



Io> ax 2 e2e "^^""du = ex ^ "^^ 



for x > xq. It is not hard to see using the inequality Dt{x) > I2 that the 
estimate from below in flT^ holds in the case where 1 + r]i < A3. 

13 



It remains to prove the estimate from below under the assumption l+r]i > 
A3. We will use the inequality Dt{x) > Ii in the proof. To estimate Ji we 
notice that z — u > z — p{z) — )■ cxd as a; — )■ 00. Therefore, Lemma [TO] can be 
applied to estimate the second integral on the right-hand side of f l30|) . This 
gives 

/p{z) 
e^Gi{u)e~'''''{z - uy^^"^ 

e^^v^exp i y^^^\ z - u)\du. 
Since the function Gi is increasing on (0, 00) and the function 

is eventually increasing, the previous inequality gives 



/i > Q;a;~2 / e2 e~^^" exp < ^ — {z — u)> du 

exp {(A3 — 1 — rji) u} du. 



Here we used the equality A^ = \ + ^|^ (see <^). Since A3 < 1 + r/i 
and p{z) — !■ 00 as 2; — !■ 00, we get Ji > aa;"^^^ ^ > xq. This establishes the 
estimate from below in Theorem [1] in the case where ^3 < 1 + 771 . 

Proof of the estimate from above in Theorem [H Let e be a small 
positive number. Denote by At{z,u) the following integral: 

y I 2ty^ 8 J 

Then formula (1271) can be rewritten as follows: 

Dt{x) = 1^^ f e^Atiz, u)pt{du) = J, + J^ + J3, (32) 

v27rtx2 Jr 

where 

Ji = . — / e^At{z,u)iJ,t{du), 

V27rtx2 J -00 
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J2 = . — / e2At{z,u)Ht{du), 

V27rtx2 Jo 



and 



J3 = r-^ — I e2At{z,u)fit{du). 

\/2TTtX2 



The number s in the previous equahties satisfies < s < 1. The value of s 
will be chosen below. 

To estimate J2, we notice that if x is large, then z — u in the expression 
for J2 is also large. Using Lemma [Eland Lemma [TOl we see that 



-^ '"'' ^" — -riiu(^_,A-'i + ^^H2V^^ 



J2 < aDt(x) + ax-^ / e^Gi(u)e''^'''(z - u)—^^^e' 



exp <; ^ (^ - u) } du. 



3 I qm 



Since the functions Gi{y) and y t-)- y^^^'^e'^^^ grow slower than the func- 
tion y H-). exp {fy} (see Lemma [9]), the previous inequality and formula ([9]) 
imply that 



sz 



J2 < ax ^^^^ + ax 2 / exp < f - - ^1 + - J "U 



< ax~'^^~^'^ + ax"'^^"^3 / exp {{A^ — 1 — %) u} du 

Jo 

< a (^— + -7^ ) (33) 

for X > xq. 

The function At is bounded (this has already been established in the 
previous part of the proof). Therefore, 

J3 < ax-i / e^Gi{u)e-'^'''du. (34) 

J sz 

Since the function Gi{u) grows slower than the function y H- exp {Cm} for 
any C > (see Lemma [9]), estimate ( l34l) implies that 



3 



J3 < ax"^+''(5+C-m)^ X > Xq. 
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Now using the fact that C can be any close to and s any close to 1, we see 
that 

•^3 < «^i:;;;rz7' ^ > ^0. (35) 

We will next estimate Ji . It follows from Lemma [8] that 

3 / u 

Ji = ax~2 / e2At{z,u)G2{u)e'^^'^du. 

J — oo 

Since m < 0, we see that 2; — m is large if x is large. Using Lemma [TDl we 
obtain 



exp l-^^^^^{z - u)\ G2{u)e'^'''du. (36) 

The function y ^^ y ^ '^e ^^^ is eventually increasing. Moreover, it grows 
slower than ez^. Since z — u> z m (l36l). we have 



Ji < ax 2 / 62 exp < ( J, V -]{z-u)\ G2{u)e'^^'^du 



00 



< ax-^^+i r et exp | f ^^^ + ^ _ f") ^| G2{u)e'^'''du 



-A3+i / I / 1 __ V8C + t , £ 



coo f 

«x-^^+ty^ exp|(-^-r/2-^^^^^ + -l«|G'2(-«)rfw. (37) 

If £ is sufficiently small, then the integral in (137|) converges (use Lemma |9]). 
It follows from (1571) that 

-^1 < tt^ — 5 X > Xq. (38) 

Finally, combining (!32|) . (133|) . (!35|) . and (138|) . we establish the estimate 
from above in Theorem [H 

Proof of Theorem [21 The following formula can be obtained from ( l22l) : 



a; / y X y V27rtx2 

e»..wf^expf-<i^iig±^-f^.,. (39, 
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It follows from (^ and (^^ that 



2; y " V a; / ^y27^ix^ 




'0 



(log ^- uf 
2ty^ 



e^Mdu) I '-^^^^ exp { - ' " -»f ^ ^ - ^ !> %, (40) 




where 



+ /" Gi(-M)e(''i-i)"dM (41) 

J Ar\{-cyofl) 



for all Borel sets A C M. In gTj), d and Ga are defined by (|25|) and (|26|), 
respectively. Now it is clear that we can use the proof of Theorem [1] with 
the pairs {rji^p) and (772,9) replaced by the pairs (772 + 1,?) and (r/i — l,j9), 
respectively. We should also take into account Lemma IHl It is not hard to 
see using fl59]) that for every e > 0, there exist constants Ci > 0, C2 > 0, and 
X > such that the following estimates hold: 

5.(ji + J^)sx-'A(^)<S.(j^ + j^) (42) 

for all X > X. The constants C2 and x depend on e. Now it is clear that ( !T6|) 
follows from ( |42|) . 

This completes the proof of Theorem [2] 

We do not include the proofs of Theorems [3] and HI because these theorems 
can be established exactly as Theorems [1] and [2j 
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